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SECTION  1 


INTRODUCTION 

Hypersonic  vehicles  during  reentry  may  experience  extensive  damage  due  to  encounters 
with  condensed  phase  H^O  in  the  form  of  either  ice  particles  or  liquid  droplets.  These  en- 
counters can  lead  to  augmented  heating  rates  and  mechanical  erosion.  The  condensed  phase  H^O 
particles,  termed  hydrorneteors , must  traverse  the  shocked  gas  enveloping  the  reentry  vehicle 
prior  to  impacting  the  surface.  The  shock  layer  can,  depending  on  flight  conditions,  vehicle 
configuration,  and  particle  size,  have  a marked  influence  on  the  hydrometeors  traversing  the 
region  between  the  shock  wave  and  the  body  surface.  Hence,  an  understanding  of  hydrometeor/ 
shock  layer  interaction  is  required  in  order  to  predict  impact  conditions  from  specified  ve- 
hicle configuration  and  free  stream  conditions. 

This  report  documents  the  results  of  an  investigation  directed  at  surveying  a portion 
of  the  phenomenology  associated  with  the  interaction  between  water  droplets  and  strong  shock 
waves.  The  investigation  has  concentrated  on  phenomena  arising  from  heat  transfer  rather  than 
from  mechanical  forces;  treatments  of  stripping  and  wave  formation  (interface  instability 
considerations)  leading  to  catastrophic  breakup  touched  upon  are  not  treated  in  depth. 

Under  many  conditions  of  interest,  all  or  a portion  of  a liquid  drop  may  attain  super- 
critical pressures  (and  temperatures)  giving  rise  to  the  question  of  what  role,  if  any,  is 
played  by  interfacial  tension  in  inhibiting  drop  breakup.  Accordingly,  the  possible  liquid 
motion  resulting  from  large  pressure  and  density  gradients  is  addressed  and  solutions  are 
presented  for  the  time  required  for  significant  droplet  distortion  under  purely  hydrodynamic 
forces. 

This  potential  breakup  mode  does  not  consider  surface  tension,  and  to  the  authors' 
knowledge  has  not  been  considered  previously.  In  Reference  1 tests  at  Mach  numbers  in  the 
1.5  to  3.5  range  on  750  to  4500  um  size  droplets,  it  was  indicated  that  the  primary  mechanism 
for  breakup  is  shear  induced  stripping.  In  Reference  2 it  is  hypothesized  that  at  higher 
Mach  numbers  catastrophic  breakup  is  dependent  on  Ueber  number. 
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For  regimes  of  interest,  the  mass  density  of  condensed  phase  materials  is  negligibly 
small  compared  to  the  mass  density  of  gas  in  the  free  stream.  Hence,  the  assumption  can  be 
made  that  the  particulate  matter  will  not  significantly  perturb  the  flow  field. 

Technology  related  to  the  present  problem  has  been  to  a certain  extent  developed  by 
investigators  concerned  with  two  phase  flow  in  combustors  and  rocket  nozzles.  References  3 
through  7 contain  related  work.  Although  the  aforementioned  references  present  results  re- 
lated to  hydrometeor  shock  layer  interactions,  the  governing  parameters  (i.e.,  particle  pro- 
perties, gas  properties,  and  flow  field  conditions)  differ  from  those  of  interest  for  this 
investigation.  Waldman  and  Reinecke  (Reference  2)  treat  the  general  problem  of  water  droplet 
shock  layer  interactions  for  reentry  flight  conditions.  The  applicability  of  the  results  in 
Reference  2 is  severely  limited  Hue  to  the  adoption  of  numerous  simplifying  assumptions.  In 
addition,  the  heat  transfer  from  the  gas  to  the  particle  is  not  treated  correctly  in  Reference 
2.  References  8 and  9 contain  corrections  to  the  formulation  of  the  heat  transfer  model  given 
in  Reference  2. 

In  References  2,  8 and  9 it  was  assumed  that  the  droplet  and  surroundings  were  in  che- 
mical equilibrium  and  that  conditions  within  the  droplet  and  the  gas  were  in  steady  state. 
Although  relaxation  times  associated  with  potential  particle/shock  interaction  phenomena  are 
small,  the  particle  residence  time  in  the  shock  layer  can  be  as  small  as  fractions  of  a micro- 
second. Hence,  it  is  necessary  to  validate  the  equi librium/steady  state  assumptions.  Once 
regimes  for  equilibrium,  steady  state  flow  have  been  defined,  mass  transfer  coefficients  can 
be  employed  to  predict  particle  mass  loss  rates.  Standard  Nusselt  number  correlations  as 
given  in  References  3 through  7 are  only  an  approximation  since  mass  transfer  effects  on  Nus- 
selt  number  were  not  included.  In  addition,  results  based  on  Nusselt  number  correlations  with 
mass  transfer  corrections  are  still  suspect  since  these  correlations  and  corrections  have  not 
been  experimentally  validated  for  the  specific  conditions  of  interest  in  the  present  investi- 
gation. Hence,  in  order  to  obtain  heat  and  mass  transfer  coefficients  with  a reasonable  degree 
of  reliability  a limited  number  of  calculations  were  performed  in  which  these  quantities  were 
obtained  from  exact  solutions  to  the  multicomponent  laminar  boundary- layer  equations. 

Previous  studies  have  indicated  that  shock  layers  will  have  negligible  effects  on  par- 
ticles in  the  rain  drop  size  range  (500  to  2500  urn);  hence,  attention  has  been  focused  on  sub- 
100  um  diameter  droplets.  The  yros^ure  and  enthalpy  ranges  considered  correspond  to  altitude 
ranges  from  20  to  40  kft,  velocity  ranges  from  10  to  20  kft/sec,  and  nosetip  radii  of  0.75  ind 
2.5  inches. 
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In  summary,  a survey  of  particle  shock  interactions  relative  to  the  following  phenomena 
have  been  made: 

• Mechanical  breakup 

• Diffusion  and  kinetic  limited  mass  loss 

• Convective  heat  and  mass  transfer 

It  is  emphasized  that  previous  studies  have  considered  the  last  item  to  be  the  significant  ther- 
mal event.  The  present  investigation  does  not  treat  the  other  phenomena  in  detail  but  focuses 
on  determining  the  flight  regimes  and  particle  dimensions  for  which  they  will  be  significant. 

The  particle  as  it  traverses  the  shock  layer  will,  at  least  initially  and  in  most  cases 
until  impact,  be  in  a doubly  shocked  gas  and  hence  an  extremely  severe  hyperthermal  environment. 
For  convenience,  a set  of  parameters  that  can  be  utilized  to  characterize  certain  pertinent 
particle  and  vehicle  parameters  has  been  generated  and  is  presented  in  Appendix  A.  These 
include  the  following  quantities  as  functions  of  vehicle  altitude  and  Mach  number: 

t Free  stream  unit  Reynolds  number 

• Sonic  point  unit  Reynolds  number 

• Vehicle  stagnation  pressure 

• Shock  standoff  (normalized  with  nose  radius) 

• Stagnation  point  momentum  thickness  (scaled  with  square  root  of  nose  radius) 

• Vehicle  total  enthalpy 

• Particle  total  enthalpy 

• Particle  Mach  number 

• Particle  Weber  number 

• Particle  stagnation  pressure 

• Particle  unit  Weber  number 
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SECTION  2 


ANALYSIS 

The  analysis  is  focused  on  characterizing  thermally  induced  events  taking  place  arising 
from  free  stream  water  droplets  interacting  with  reentry  vehicle  shock  layers.  As  was  men- 
tioned, mechanical  breakup  arising  from  high  induced  internal  pressures  is  also  considered. 

The  particle  size  range  considered  in  these  analyses  is  smaller  than  that  generally  encountered 
in  precipitating  rain  but  is  representative  of  that  found  in  cloud  formations.  It  is  empha- 
sized that  the  particle  environment  consists  of  a doubly  shocked  region  and  hence,  in  terms  of 
enthalpy  and  pressure,  it  will  be  more  severe  than  the  vehicle  shock  layer.  The  range  of  con- 
ditions employed  in  the  analyses  is  summarized  below. 

Configuration  Condition 

Nose  radii:  0.75  and  2.50  in. 

Free  Stream  Conditions 
Altitude:  15-40  kft 
Velocity:  10-20  kft/sec 

Particle  Conditions 

Rp:  5-100  urn 

P$:  50-2000  atm 

Ht:  5000  - 24,000  Btu/lbm 

M : 2.4  - 2.9 

P 

Appendix  A gives  further  details  relating  free  stream  conditions  to  particle  and  nosetip  en- 
vironmental conditions. 

In  Section  2.1,  the  entry  of  particles  in  a shock  layer  is  described  and  a simple  mo- 
del is  used  to  compute  droplet  distortions.  Section  2.2  compares  the  times  to  achieve  equi- 
librium mass  diffusion  with  particle  residence  times  and  shows  that  convective  film  coefficient 
modeling  in  the  shock  layer  is  appropriate.  Section  2.3  presents  analyses  of  the  convective 


heat  and  mass  transfer  between  a particle  and  the  doubly-shocked,  shock  layer  gases  including 
particle  transient  thermal  response  calculations. 

2.1  MECHANICAL  BREAKUP 

Mechanical  forces  are  present  in  the  shock  layer  of  the  vehicle  (or  induced  in  it) 
which  tend  to  cause  asymmetric  loads  on  the  water  droplet  and  can  cause  it  to  break  up.  In- 
itially, the  droplet  encounters  asymmetric  loads  as  it  crosses  the  shock  front.  Subsequently, 
a quasi -steady  gasdynamic  flow  is  established  about  the  droplet  introducing  shear  forces  along 
its  surface  (which  are  viewed  as  a stripping  effect  and  not  considered  here)  and  a pressure 
difference  between  the  stagnation  point  and  other  points  on  the  surface  of  the  droplet.  The 
droplet  cannot  support  such  asymmetric  loads  and  must  respond  hydrodynamically.  If  sufficient 
time  is  available  before  the  droplet  impacts  the  surface  of  the  body,  it  can  distort  to  the 
point  where  it  can  be  viewed  as  having  broken  up.  The  sequence  of  events  is  shown  in  Figure 
2-1,  which  were  inferred  from  theoretical  considerations. 


Figure  2-1.  Droplet  breakup  phenomenology. 


Entry  Through  the  Shock  Front 

As  the  droplet  enters  the  shock  laye*-  jf  the  body,  it  displaces  the  shock  front  lo- 
cally and  simultaneously  initiates  a new  gasdynamic  shock  wave  inside  the  shock  layer  and  a 
compression  wave  inside  the  droplet  itself  (Figures  2-la  and  2-lb).  As  the  speed  of  sound 
in  water  is  typically  about  5000  ft/sec,  or  1/3  to  1/2  the  velocity  of  the  droplet  relative 
to  the  shock  layer  ga^es,  the  compression  wave  in  the  water  will  not  have  time  to  make  a full 
traverse  across  tb*  droplet  before  the  latter  is  completely  through  the  shock  front  and  into 
the  shock  layer.  It  follows  that  there  will  be  a region  within  the  droplet  where  there  is 
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undisturbed  water.  Thus,  water  droplets  will  always  be  able  to  pass  through  the  shock  front 
without  breaking  up,  provided  only  that  the  droplet  velocity  is  greater  than  the  sound  speed 
of  water. 


Traverse  of  the  Shock  Layer 

A quasi-steady  gasdynamic  flow  will  be  established  about  the  droplet  very  quickly  after 
the  latter  enters  the  vehicle  shock  layer.  This  supersonic  flow  will  induce  significant  hy- 
drodynamic pressure  differences  across  the  droplet.  Of  course,  such  pressure  differences  re- 
present a potential  driving  force  for  flow  within  the  droplet.  Internal  motion  would  cause 
distortions  and  possibly  break  up  as  indicated  in  Figures  2-ld  and  2-Ie.  The  magnitude  of 
this  effect  can  be  estimated  by  the  use  of  simple  force  laws. 

Consider  a droplet  of  spherical  shape  of  initial  radius  Rp,  density  p,  stagnation  pres- 
sure P$  and  a pressure  on  its  back  side  which  can  be  neglected  relative  to  P$.  The  net 

force  tending  to  accelerate  the  particle  is  then 

F-P$nR’  (2-1) 


while  the  mass  to  be  accelerated  is 


M 


P 


(2-2) 


If  the  acceleration  acts  over  a distance  ARp,  which  can  also  be  viewed  as  the  increase  in 
radius  of  the  particle  due  to  distortion,  then 


4Rp 


2 

res 


(2-3) 


where  a is  the  acceleration  and  t is  the  residence  time  of  the  droplet  in  the  shock  layer 

res 

of  the  body.  Then  combining  the  abo'e  equations  with 


a 


(2-4) 


yields 


p P2 

MP  P 


res 


(2-5) 
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which  is  the  key  result  desired.  The  residence  time,  t^,  is  fixed  by  the  nose  radius  of  the 
vehicle  and,  consequently,  does  not  vary  significantly.  Likewise,  the  water  density  does  not 
vary  significantly,  so  that  breakup  is  primarily  a function  of  droplet  radius  and  stagnation 
point  pressure.  Indeed,  distortion  (and  hence,  the  likelihood  of  breakup)  is  seen  to  increase 
with  the  stagnation  pressure  and  decrease  with  the  size  of  the  initial  droplet  radius,  in 
rough  agreement  with  t1  available  experimental  data.  It  is  viewed  as  significant  that  sur- 
face tension  does  not  play  an  important  role  in  this  model,  in  contrast  to  tne  model  of  Wald- 
man  and  Reinecke  (Reference  2) 

A number  of  solutions  to  Equation  (2-5)  are  presented  in  Table  2-1,  considering  in- 
teresting flight  conditions  and  for  a residence  time  of  2 x 10"7  seconds.  As  shown  there, 
the  10  urn  droplets  are  likely  to  break  up  for  all  the  interesting  pressures,  as  are  the  20  pm 
droplets.  The  larger  droplets  will  probably  not  break  up,  even  at  the  highest  pressures  of 
interest. 


TABLE  2-1.  DROPLET  DISTORTIONS 


Pressure  (atm) 

500 

1000 

2000 

Radius  (urn) 

£Rp/Rp 

10 

3.75 

7.5 

15. 

20 

0.94 

1.88 

3.76 

50 

0.149 

0.298 

0.596 

100 

0.0371 

0.0752 

0.1504 

2.2  DIFFUSION  AND  KINETIC  LIMITED  MASS  LOSS 

If  it  is  assumed  that  the  hydrometeor  surface  is  in  vaporization  equilibrium  and  radi- 
ation can  be  neglected,  an  established  correlation  exists  between  a parameter  representing 
nondimensional  blowing,  B^,  and  an  appropriate  combination  of  total  enthalpy,  He»  surface  en- 
thalpy, Hw,  and  heat  of  formation  of  surface  material,  Hf.  The  correlation  is  given  by* 


6'  ■ f.n 
o 


' *> 
V*- 


-6) 


*This  relation  is  identical  to  relations  employed  in  Appendix  C. 
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where 


b; 


"d 


VST 


(2-7) 


The  quantity  ru  is  the  diffusion  limited  surface  mass  flux  normal  to  the  surface  and  p u is 
u e e 

the  mass  flux  parallel  to  the  surface  in  the  inviscid  flow.  In  order  to  assess  the  signifi- 
cance of  transient  effects  in  limiting  diffusion  mass  transfer,  the  time  constant  for  this 
phenomenon  must  be  compared  to  the  particle  residence  time  in  the  shock  layer.  The  time  con- 
stant appropriate  for  diffusion  through  the  viscous  layer  is  given  by  the  ratio  of  mass  in  the 
boundary  layer  to  diffusion  mass 


(2-8) 


The  integral  lor  the  boundary  layer  mass  in  the  previous  equation  can  be  approximated  by  (Re- 
ference 10) 


(2-9) 


and  the  value  mg  is  obtainable  from  Equations  (2-6)  and  (2-7)  for  specified  flight  conditions. 

Over  the  conditions  of  interest  it  can  be  shown  that  the  quantity  (p^C^)  Is  within 
the  range  200  to  1900  lbm/ft2sec.  The  quantity  tQ  is  calculated  for  the  high  end  of  the 
droplet  pressure  range,  which  will  be  on  the  order  of  2000  atm  (see  Appendix  A).  Table  2-2 
contains  time  constants  for  the  range  of  conditions  of  interest. 

TABLE  2-2.  DIFFUSION  TIME  CONSTANT  tQ  (SEC  x 10’)  AT  P * 2000  ATM 


(Btuj^bm) 

peueCHo  ^bm/f^sec) 

200 

500 

1000 

5,000 

40 

7 

1.7 

10,000 

20 

3 

0.8 

15,000 

8 

1.3 

0.3 
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Particle  shock  residence  times  for  the  stagnation  region  are  of  the  order  of  2 x 10‘7 
seconds,  hence,  from  Table  2-2  It  does  not  appear  that  transients  are  Important  In  determin- 
ing diffusion  limited  mss  loss. 

Once  the  steady  boundary  layer  becomes  filled  with  vaporization  products,  mass  loss 
will  be  diffusion,  rather  than  kinetic  controlled.  An  estimate  of  the  times  associated  with 
vaporization  kinetic  events  can  be  made  based  on  simple  considerations.  A rate  equation  of 
the  form 


\ * f(V  - kr  Pw  <2-10> 

can  be  written  to  describe  the  heterogeneous  vaporization  process.  The  rate  coefficient,  k^, 
can  be  evaluated  from  simple  kinetic  theory  as  follows: 

kr  * 122  x 106  lbm/ft2sec-atm  (2-11) 

This  expression  includes  a sticking  factor  of  units  which  has  an  uncertainty  of  perhaps  a 
factor  of  25.  Using  the  steam  table  (Reference  11)  and  detailed  balancing,  the  net  forward 
rates  can  be  estimated  as  follows: 


Temperature,  T* 
(°R) 

Vapor  Pressure 
(atm) 

(lb/ft^sec) 

575. 

0.1 

2 x 10‘ 

672. 

1.0 

2 x 107 

817. 

10.0 

2 x 10* 

Pressures  of  actual  Interest  are  much  higher  than  given  above,  but  the  trend  to  higher  m^ 
values  with  increasing  pressure  Is  clear.  Above  the  critical  pressures,  vaporization  kinetics 
modeling  Is  not  clear  but  since  the  process  Is  homogeneous,  faster  rates  should  occur. 

The  time  required  to  supply  the  boundary  layer  with  vaporization  products  Is  a measure 
of  the  time  when  vaporization  kinetics  would  be  signfl leant.  This  time,  t^,  may  be  estimated 
from 


6 

*r  * </  °bid>:/\ 
0 


(2-12) 
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Typical  values  for  heat  transfer  coefficient  values  of  Interest  are  given  In  Table  2-3. 
TABLE  2-3.  KINETIC  TIME  C INSTANT  tK  (SEC  xlOn)tP«  2000  ATM 


Ts 

<°R) 

(ififtfcc) 

200 

500 

1000 

575 

50 

20 

10 

672 

5 

2 

1 

812 

0.5 

0.2 

0.1 

Again  since  hydrometeor  residence  times  are  on  the  order  2 x 10"’  seconds,  it  Is  safe  to  as- 
sume that  vaporization  kinetic*  are  not  significant. 

The  above  analyses  mean  that  continuum,  equilibrium,  fluid  mechanics  apply  in  the 
vir.inity  of  a hydrometeor  and  that  boundary  layer  modeling  can  be  used  to  investigate  hydrometeor/ 
shock  layer  Interactions.  Film  coefficient  modeling  of  heat  and  mass  transfer  events  is  described 
is  the  following  section. 

2.3  CONVECTIVE  HEAT  AND  MASS  TRANSFER 

The  air  flowing  about  the  surface  of  a water  droplet  (Immersed  in  a shock  layer)  is  sub- 
jected to  shock  heating  by  two  successive  near-normal  shock  fronts.  The  levels  of  pressure  and 
enthalpy  which  result  give  the  air  the  potential  for  delivering  exceedingly  high  convective 
heating  rates  to  the  droplet  surface.  If  such  heating  rates  are,  in  fact,  delivered,  one  would 
expect  significant  mass  removal  from  the  droplet  by  ablation.  The  prediction  of  the  coupled 
convective  heating-ablative  response  is  made  complicated  and  uncertain  by  a number  of  factors 
not  usually  encountered  in  conventional  problems. 

The  novel  features  of  this  problem  stem  from  the  fact  that  the  water  can  behave  much 
like  a gas  at  the  supercritical  pressures  of  Interest.  Therefore,  a free  shear  layer  will 
exist  along  the  air-water  boundary  in  which  both  the  air  and  water  co-exist  and  are  responding 
mechanically.  Questions  immediately  arise  with  regard  to  the  proper  treatment  of  the  physical 
events.  It  is  Important  to  employ  a solution  procedure  which  couples  the  motion  of  the  water 
to  that  of  the  air?  Alternately,  can  the  droplet  surface  be  approximated  as  being  stationary? 

What  are  the  parameters  which  govern  the  problem?  Further,  as  there  is  no  phase  change  between 
the  water  in  the  droplet  and  the  water  which  has  diffused  into  the  surrounding  air,  how  does 
one  decide  if  mass  ha;  been  removed  from  the  droplet  or  not?  Also,  one  would  expect  to  have 
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a significant  amount  of  mass  removed  by  stripping  if  existing  theories  are  ;o  be  believed  (Re* 
ference  2).  Will  mass  transfer  effects  significantly  reduce  the  stripping  effect?  Indeed,  is 
the  mass  stripped  from  the  front  half  of  the  droplet  removed  from  the  droplet  or  merely  moved 
to  its  backward- facing  half?  Finally,  when  is  it  appropriate  to  consider  steady-state  fluid- 
dynamic  and  ablative  response  events,  and  when  must  transient  events  be  considered? 

In  the  material  presented  in  the  following  paragraphs  a number  of  the  questions  raised 
with  regard  to  proper  modeling  approaches  are  discussed.  The  objective  of  this  material  is 
to  reduce  this  exceedingly  complex  problem  to  a point  where  valid  estimates  of  the  ablative 
mass  losses  can  be  obtained. 

Inviscid  Flow  Parameters 

The  inviscid,  adiabatic  flow  field  events  of  Interest  include  the  pressure  distribu- 
tion and  boundary  layer  edge  velocity  and  enthalpy.  The  results  of  an  "exact"  solution  were 
used  to  obtain  the  pressure  distributions  assuming  an  invariant  spherical  shape.*  The  boun- 
dary layer  edge  velocities  were  obtained  by  calculating  an  isentropic  expansion  to  the  local 
pressure.  Finally,  the  stagnation  point  pressures  and  enthalpies  were  obtained  using  a pro- 
cedure (the  ACE  program.  Reference  12)  which  solves  the  Rankine-Hugoniot  equations  governing 
the  double  shock  system,  coupled  to  the  equations  governing  the  equilibrium  chemical  state. 

The  free  stream,  downstream  (after  the  Initial  shock),  and  stagnation  point  solutions  for  the 
flight  conditions  of  interest  are  presented  in  Tables  2-4,  2-5,  and  2-6,  respectively.  Addi- 
tional material  on  the  inviscid  flow  will  be  presented  and  discussed  with  the  discussion  of 
the  convective  heating  events. 

Air/Water  Mixing  Region 

The  air/water  mixing  region  ultimately  determines  the  amount  of  mass  removed  from  the 
front  surface  of  the  droplet;  consequently,  a solution  to  its  governing  equations**  would  re- 
present a major  step  toward  the  solution  of  the  ablation  problem.  To  the  authors'  knowledge 
solutions  of  this  type  are  not  available  in  the  literature.  Moreover,  the  generation  of  such 


Of  course,  this  represents  an  approximation  as  the  droplet  shape  will  change  (in  general) 
as  it  traverses  the  shock  layer.  Shape  change  effects  could  have  been  considered  but  it 
was  felt  that  they  would  greatly  complicate  tne  analysis  without  adding  appreciably  to  the 
understanding  of  the  essential  physical  events. 

**The  mixing  region  is  governed  by  the  boundary  layer  equations  (global  mass,  momentum,  energy 
and  diffusion)  which  must  be  solved  subject  to  both  a valid  set  of  transport  properties  and 
to  a chemical  equilibrium  constraint  valid  'rom  both  the  pure  gas  phase  (air)  into  the  pure 
liquid  phase  (water). 
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TABLE  2-4.  FREESTREAM  PROPERTIES 


) (Btu/lb)  where  Tref  = 536.698°R  (298.16°K) 


TABLE  2-5.  DOWNSTREAM  PROPERTIES 
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TABLE  2-6.  PARTICLE  SHOCK  LAYER  PROPERTIES 


where 


solutions  is  beyond  the  scope  of  the  present  study.  Fortunately,  problems  have  been  studied 
(and  solutions  published)  which  include  many  of  the  essential  features  of  the  present  problem 
Solutions  to  such  problems  give  qualitative  features  and  are  useful  in  guiding  the  construc- 
tion of  approximate  solution  procedures. 

R.  C.  Lock  (Reference  13)  considered  the  two-dimensional  motion  of  a stream  of  fluid 
with  velocity,  U1 , density,  pj,  and  viscosity,  u1 , over  a parallel  stream  with  velocity,  U2* 
density,  p2»  and  viscosity,  u2*  Both  fluids  were  assumed  to  be  incompressible.  The  configu- 
ration is  shown  in  Figure  2-2. 


The  equations  and  boundary  conditions  are  written  out  in  Appendix  B.  Self-similar  solutions 
were  obtained  in  terms  of  the  coordinate 


n 


1 


(2-13) 


for  the  upper  stream  and  the  coordinate 


(2-14) 


for  the  lower  stream. 

Numerical  solutions  to  Lock’s  problem  are  presented  in  Figure  2-3  for  the  interesting 
case  of  U2  = 0 and  for  a range  of  viscous  density  ratios  [z  p^/p^).  It  is  significant 
that  the  magnitude  of  the  viscous  density  ratio  alone  determines  the  importance  of  the  coupling 
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Figure  2-3.  Solutions  to  Lock's  problem. 
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between  the  two  streams.  Thus,  for  values  of  this  parameter  greater  than  or  equal  to  roughly 
100,  one  should  be  able  to  approximate  the  upper  stream  as  a boundary  layer  over  a stationary 
wall.  Indeed,  if  mass  transfei  affects  were  included,  this  number  would  probably  be  somewhat 
lower.  This  result  should  also  be  applicable  to  the  air/water  mixing  region  over  the  droplet; 
that  is,  if 


pwater1Wer/‘WJair  ' 0(100) 

the  motion  of  the  water  can  be  uncoupled  from  the  motion  of  the  air. 

Typical  values  of  density  and  viscosity  are  presented  in  Table  2-7  for  air  edge  condi- 
tions and  for  water.  As  shown,  the  viscous  density  ratio  is  often  -0(100),  thereby  justify- 
ing uncoupling  the  water  flow  from  the  air  flow.  Based  on  this  conclusion,  unblown  convective 
heating  rates  were  computed  from  boundary  layer  theory  for  conditions  and  droplet  sizes  of  in- 
terest. These  are  discussed  below. 


TABLE  2-7.  DENSITIES  AND  VISCOSITIES 


Air 

P (atm) 

p,  (lb/ft3) 

Uj  (lb/ft-sec) 

PI  PI 

943 

1.197 

1.244  x 10"1* 

1.5  x 10'4 

1900 

2.155 

1.313  x 10"4 

2.83  x 10  14 

2910 

3.027 

1.366  x 10-14 

4.1  x 10"4 

Water 

T (°K) 

P2  (lb/ft3) 

1*2  (lb/ft-sec) 

J2^2_ 

50 

62.4 

3.66  x 10"14 

228.  x lO-4 

100 

62.4 

1.9  x 10”4 

117.  x 10‘4 

Unblown  Convective  Heating  Rates 

The  analysis  was  limited  to  spherical  droplets  traversing  the  shock  layer  on  the  vehicle 
stagnation  streamline.  Solutions  were  obtained  to  the  laminar,  nonsimilar,  chemically  reacting 
boundary  layer  equations  using  a detailed,  numerical  procedure  (the  BLIMP  code.  Reference  14). 
Vehicle  free  stream  conditions  encompassed  one  altitude  range  from  15-40  kft  and  vehicle  Mach 
numbers  from  10  to  20.  The  vehicle  free  stream  conditions  were  taken  from  Reference  15  and  are 
summarized  in  Table  2-4.  Conditions  behind  the  vehicle  shock  were  obtained  from  normal  shock 
ACE  solutions  (Reference  12),  and  are  summarized  in  Table  2-5.  The  free  stream  conditions 


behind  the  particle  shock  were  also  obtained  from  normal  shock  ACE  solutions  using  conditions 
behind  the  vehicle  shock  as  the  properties  upstream  of  the  particle  shock.  These  results, 
which  are  used  as  the  free  stream  conditions  for  the  BLIMP  calculations,  are  summarized  in 
Table  2-6.  Note  that,  for  the  flight  conditions  considered,  particle  Mach  numbers  (M^)  lie  in 
the  narrow  range  from  2.4  to  2.9. 

Pressure  distribution  around  the  sphere  was  obtained  from  an  "exact"  solution  to  the 
inviscid  flow  field  around  a sphere.  The  ratio  of  local  static  pressure  to  stagnation  pres- 
sure (P/P$)  is  shown  in  Figure  2-4  for  an  upstream  Mach  number  of  1.8.  The  flow  field  over 
the  rear  hemisphere  is  not  well -characterized  since  separation  occurs  in  the  region  90°  < 0 < 
120°.  Therefore,  heating  predictions  were  only  made  over  the  windware  side  of  the  sphere, 
i.e. , 0°  < 0 < 90°. 

Heat  transfer  distributions  were  calculated  for  droplet  sizes  from  5 to  100  ym  at  15 
kft  and  40  kft  for  free  stream  Mach  numbers  of  10,  15,  and  20.  Predicted  distributions  of 
cold-wall  heat  transfer  coefficient  are  shown  in  Figure  2-5  for  altitudes  of  40  kft  and  15 
kft,  respectively.  These  results  were  correlated  and  the  similarity  of  the  distributions 
suggested  a curve  fit.  The  heat  transfer  distribution  can  be  represented  in  terms  of  the 
stagnation  point  values  as: 

peueCHQ  j cos  6 for  0°  < 0 < 60° 

(pelleCH0)$P  l j exp[-( 0-601/(90-60)]  for  60“  < 0 < 90° 

The  stagnation  point  heat  transfer  was  correlated  and  can  be  represented  by: 

[peueCH  3SP  = 33^M«’-” 

0 

where  [p  u C[\  ]CD  ■ stagnation  point  heat  transfer  coefficient,  lbm/ft2sec 

6 6 0 Dr 

P = free  stream  static  pressure,  atm 
00 

Rp  = droplet  radius,  um 

The  applicability  of  the  correlation  was  shown  by  comparing  the  nondimensional  distribution 
to  BLIMP  results  for  a droplet  with  a radius  of  50  ym  at  25  kft  for  Mach  numbers  of  15  and 
20.  The  results  of  this  comparison  is  shown  in  Figure  2-6.  Estimates  of  the  heat  transfer 
In  the  separated  region  range  from  1 to  4 percent  of  stagnation.  The  average  heat  transfer 
over  the  front  half  of  the  sphere  is  71  percent  of  stagnation. 
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Figure  2-4.  Sphere  pressure  distribution. 


2-6.  Comparison  of  heat  transfer  distribution  correlation  result  to  BLIMP 
result  for  altitude  = 25  kft,  particle  radius  * 50  ym. 


Values  of  the  cold  wall  heating  rates  are  also  useful  and  are  tabulated  along  with 
transfer  coefficient  and  total  enthalpy  values  in  Table  2-8  for  the  droplet  stagnation  point 
conditions.  Note  that  the  rates  are  in  the  range  2 x 105  to  10  x 10*  Btu/ft2sec  which  are 
huge  relative  to  those  encountered  by  reentry  vehicles  in  even  the  most  severe  reentry  environ- 
ments . 


Transient  Heat  Conduction 

Despite  the  significant  heating  rates  experienced  by  the  particle,  vaporization  will 
not  take  place  if  the  surface  temperature  is  below  its  saturation  value.  It  is  pointed  out 
that  conditions  including  pressure,  and  hence  saturation  temperatures,  vary  around  the  droplet 
during  its  traverse  of  the  shock  layer.  Moreover,  prediction  of  this  variation,  especially 
past  separation,  is  exceedingly  difficult  and  matters  are  complicated  by  the  existence  of  a 
supercritical  state  for  certain  flow  conditions.  Since  the  present  calculations  are  designed 
to  determine  for  which  flight  regimes  various  mechanisms  need  be  considered,  approximations 
are  made.  Thus,  the  time  to  attain  the  saturation  temperature  for  P * P$,  used  as  a measure 
of  the  surface  relaxation  time  for  vaporization.  It  is  also  assumed  that  the  hydrometeors  are 
spherical. 

The  unsteady  heat  conduction  solution  for  the  surface  of  a heated  sphere  initially  at 
T^  exposed  to  a temperature  T£  is  given  by 


where 


y t>  - t, 

T2-T« 


- = f(Fo,Bi) 

(2-17) 

ot/S' 
w p 

(2-18) 

hR 

'T* 

(2-19) 

The  heat  transfer  coefficient,  h,  can  be  obtained  from  the  stagnation  point  convective  heating 
correlation  given  as  Equation  (2-16).  That  is. 


CP,air  peueCH 


p.air 


(2-20) 
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TABLE  2-3.  STAGNATION  POINT  COLD  MALL  HEATING  RATES  {qQ) 


z 

(kft) 

M 

oo 

RP 

(ym) 

% 

(Btu/ft2sec) 

[peueCHo^5P 

(lb/ft2sec) 

Htot 

(Btu/lbm) 

4 

■ 

a 

5 

4,431,000 

327 

13,563 

■ 

20 

2.216,000 

163 

1 

1 

100 

990,900 

73 

20 

113 

7,516 

D 

5 

129 

3,228 

25 

l 

50 

2,207,000 

147 

14,964 

| 

15 

1 

830,600 

104 

8,459 

15 

2 

a 

5 

9,412,000 

580 

16,230 

20 

4,706,000 

290 

50 

2,977,000 

183 

100 

2,105,000 

130 

1 

5 

5 

3,583,000 

399 

8, 

983 

20 

1,791,000 

199 

, 

m 

5 

903,400 

232 

3,898 

U 

100 

202,000 

51 

i 
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Heat  transfer  coefficient  values  given  by  the  above  analysis  Indicate  Nusselt  numbers  which 
are  2 to  4 times  the  values  given  by  typical  low  Mach  number  correlations  such  as  those  used 
In  Reference  2,  8,  and  9.  If  the  Fourier  number  (Fo)  Is  very  small,  the  solution  represented 
by  Equation  (2-14)  assumes  the  semi-infinite  plate  form.  That  is,  for 

Fo  « 1 


T.(t)  - T 


T2-T. 


g(Bi  /Fo) 


(2-21) 


It  is  pointed  out  that  residence  times  In  the  shock  layer  will  be  of  the  order  of  a microsec- 
ond; the  thermal  diffusivity  a of  water  can  be  taken  as  1.8  x 10“6  ft2/sec;  hence,  for  drop- 
lets on  the  order  of  10  ym,  the  value  of  Fo  is  of  th«>  order  10“3.  For  flight  conditions  of 
Interest  the  Biot  numbers  (Bi)  fall  generally  within  th*»  following  range: 

1 i Bi  < 100 


Values  of  Fourier  number  will  be  not  greater  than  10*2.  The  functions  f(Bi,Fo)  and  g(B1,/Fo) 
satisfying  Equations  (2-17)  and  (2-21)  are 


n=l 


sin  a„  “ cos  a 
n n 


udn  * an  cos  an  s1n 


n 

a„_ 


-a2 

(sin  an)e  n 


where  a„  are  roots  of  a cot  a - (Bi  - 1) 
n n n 


(2-22) 


9 = 1 - [1  - erf(B1,^Fc)]eBl2Fo  (2-23) 

Comparisons  of  values  for  the  exact  solution  given  by  Equation  (2-22)  to  the  semi-infinite 
slab  solution  given  by  Equation  (2-23)  are  shown  in  Table  2-9. 


TABLE  2-9.  VALUES  OF  (T$  - TJ/(T‘  - T J 


Bi 

Fo 

Exact 

Solution 

Semi-Infinite 
Plate  Soi‘n 

% Error 

1 

0.010 

0.1128 

0.1035 

8.24 

1 

0.001 

0.0357 

0.0347 

2.70 

100 

0.010 

0.9528 

0.9414 

1.19 

100 

0.001 

0.3362 

0.8292 

0.83 
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These  results  indicate  the  adequacy  of  the  simpler  approximate  solutions;  thus  Equation  (2-23) 
is  employed  to  yield  the  surface  temperature  time  relationship  for  droplets  in  the  shock  layer. 
For  most  conditions  of  interest  this  equation  will  be  accurate  to  within  less  than  1 percent. 

It  is  pointed  out  that  the  heat  transfer  coefficient  correlation,  given  previously, 
requires  no  mass  transfer  correction  since  the  present  calculation  is  to  establish  times  re- 
quired to  reach  vaporization  temperatures;  thus,  this  portion  of  the  analysis  is  related  to 
heat  transfer  effects  prior  to  substantial  vaporization. 

Surface  temperature  T is  computed  using  T s T at  t * 0 as  an  initial  condition, 
s s * 

Values  are  to  be  compared  to  time  intervals  corresponding  to  the  particle  residence  time,  t , 
in  the  shock  layer,  which  is  given  by 

tres  * gV”  CexP(SAs)  " (2-24) 

00 


where 


In  the  analysis,  a hypersonic  approximation  for  the  shock  standoff  distance  is  employed,  so 
that 


As/Rn  - 0.08 


Also 


Based  on  these  assumptions,  the  approximate  residence  tirae<  of  the  hydrometeor  in  the  shock 
layer  along  the  vehicle  stagnation  line  are  as  follows 


R 

R 


N 

N 


Residence  Times  - 



- usee 

M - 10 

TO 

IS 

H.  - 20 

0.6 

0.4 

0.3 

1.9 

1.3 

1.0 
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Variations  In  altitude  and  drop  size  over  the  range  of  Interest  changes  the  above  values  less 
than  +20  percent. 

The  time  required  for  a particle  to  achieve  a condition  where  vaporization  commences 
must  be  less  than  the  residence  time  for  convective  mass  loss  to  be  significant.  Since  the 
droplets  are  typical ly  at  supercritical  pressures  for  flight  condition  of  Interest,  the 
vaporization  condition  Is  not  well  defined,  for  the  purpose  of  this  simple  analysis,  vapori- 
zation Is  assumed  to  commence  when  droplet  surface  temperature  reaches  the  critical  point  tem- 
perature, T$  * U65°R.  For  a particular  altitude  and  Mach  number  condition  the  left-hand 
side  of  Equation  (2-21)  can  be  evaluated.  Then,  Equations  (2-18),  (2-19),  and  (2-20)  can  be 
used  to  evaluate  the  time,  t,  to  achieve  vaporization. 

Figure  2-7  presents  the  times  to  achieve  vaporization  conditions  as  a function  of  the 
droplet  radius  for  a wide  range  of  altitude  and  Mach  number  conditions.  Figure  2-8  shows  the 
time  to  achieve  vaporization  as  a function  of  local  heat  transfer  coefficient.  Transfer  coef- 
ficient partially  correlates  the  results.  A further  collapse  Is  obtained  by  correlating  with 
the  stagnation  point  heat  flux.  That  Is 


DeueCH0'sP 


M2 


(2-25) 


Based  on  this  approximation,  time  to  reach  critical  temperature  at  the  droplet  stagnation 
point  is  given  by 


1050  R 

tCr1t  ' 1050  ^ wee  (2-26) 

P M 
00  00 

where  Pw  * free  stream  pressure,  atm 

Rp  • Initial  droplet  radius,  urn 

M ■ free  stream  Mach  number 
00 

t . » has  units  of  usee 
crlt 

Clearly  the  critical  time  Is  quite  sensitive  to  the  vehicle  Mach  number.  As  Figures  2-7  and 
2-8  show,  critical  times  to  reach  vaporization  conditions  at  the  droplet  stagnation  point  are 
well  below  droplet  residence  times  for  typical  reentry  flight  conditions 
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Droplet  radius  - urn 


Figure  2-7 


Droplet  vaporization  initiation 


2-25 


temperature 


It  should  be  pointed  out  that  previous  experimental  work  by  AVCO  (Reference  16)  and 
Aerotherm  (Reference  17)  was  for  snock  velocities  not  exceeding  11.7  kft/sec.  This  work  in- 
dicates that  the  principle  shock/droplet  interaction  phenomena  consists  of  droplet  accelera- 
tion, stripping,  and  breakup.  Inability  to  detect  vaporization  in  these  experiments  of  100+ 
urn  particles  can,  according  to  the  present  analysis,  be  partially  attributed  to  the  fact  that 
surface  temperatures  did  not  reach  vaporization  levels.  Consequently,  although  the  aforemen- 
tioned experiments  do  not  indicate  that  vaporization  is  taking  place  it  cannot  be  concluded 
that  the  effect  of  this  phenomenon  is  unimportant  for  all  flight  conditions. 

Mass  Transfer  Effects  on  Convective  Heating 

In  the  event  that  conditions  are  such  as  to  produce  vaporization,  the  effects  of  blow- 
ing on  convective  heating  should  be  considered  in  the  evaluation  of  mass  loss.  The  neat 
blockage  effects  have  traditionally  correlated  through  a blowing  correction  to  the  Stanton 
number  of  the  form  (Reference  13). 


ch  , 2ab; 

CH0  exp(2*B^)  - 1 


(2-27) 


where 


B 


o 


peueCH, 


(2-28) 


and  where  A is  a parameter  with  a value  that  is  a function  of  whether  the  flow  is  laminar  or 
turbulent  and  is  a function  of  the  molecular  weight  of  the  ablated  vapors.  The  results  of  a 
study  by  Bartlett,  Nicolet  and  Howe  (Reference  19)  are  presented  in  Figure  2-9,  where  it  is 
shown  that  the  correlation  is  excellent,  the  effect  of  molecular  weight  is  second  order,  and 
a value  A = 0.6  does  well  in  correlating  the  laminar  data.  At  significantly  higher  blowing 
rate.,  which  are  characteristic  of  water  droplet  vaporization,  a higher  order  of  A would  be 
applicable.  However,  since  a higher  value  would  minimize  droplet  mass  loss,  the  0.6  value  is 
used  in  subsequent  numerical  work  to  estimate  the  upper  limit  of  droplet  mass  loss. 

Droplet  Mass  Loss  by  Ablation 

The  coupled  drag  and  ablation  of  a droplet  was  computed  for  the  conditions  of  interest. 
The  approach  of  Jaffe  (References  8 and  9)  was  modified  to  include  the  effects  of  blowing  and 
real  gas  properties  on  the  heat  and  mass  transfer  about  the  droplet.  Appendix  C presents  the 
derivation  and  discussion  of  the  hydrometeor  mass  less  and  velocity  equations  which  account 
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number  due  to  surface  mass  addition  (from  Reference  11) 


for  the  above  effects.  The  relations  for  droplet  velocity  and  mass  after  the  droplet  passes 
through  the  shock  layer  are  as  follows: 


where 


Uf  = exp 


2Kn 


(1  - \V1  - 3/2  K') 

L 4 E. 


} 0 < K'  < 2/3 


(2-29) 


M„  . = 11  (1  - 3/2  KI)2 

p,f  p,i  E 


KD  = 8 C0 


(pj\  A. 

VP/  RP.i 


4' 


°-71  [peueCH0]SPh  *"  ^ + 2>  "tot'"? 

K - ¥ pp^ 


In  the  above  relations  subscript  i denotes  conditions  existing  on  the  particle  at  initial  en- 
try in  the  shock  layer.  For  example,  [PeuecH03$p  -j  indicates  the  initial  particle  stagnation- 
point  heat  transfer  coefficient  as  it  enters  the  shock  but  prior  to  any  acceleration  or  abla- 
tion occurs.  The  transfer  coefficient  can  be  obtained  from  Equation  (2-16)  for  R = R . . 

P P 

ic  that  K£  = 2/3  corresponds  to  complete  ablation  before  the  droplet  reaches  the  nosetip, 

;d  !'£  = 2/3  determines  the  critical  conditions  for  no  hydrometeor  impacts. 

Assumptions  contained  within  Equations  (2-29)  are  as  follows: 

• Stagnation  point  heat  transfer  scales  with  u//R^  as  derived  by  Bartlett  and  Putz 
(Reference  20) 

t Heat  transfer  to  the  back  of  the  droplet  is  negligible 

• Average  heat  flux  to  the  front  of  the  droplet  is  71  percent  of  the  stagnation  point 
flux 

• Rate  of  energy  tenge  in  the  droplet  is  negligible  cor  pared  to  convective  and 
vaporization  energy  rates. 

Other  details  of  the  formulation  are  given  in  Appendix  C. 

Figure  2-10  shows  hydrometeor  mass  reduction  results  for  l-lach  20  conditions  as  a func- 
tion of  shock  standoff  distance.  For  nose  radii  and  drop  dimensions  of  interest  in  this  study, 
final  mass  is  always  more  than  85  percent  of  initial  mass  even  for  very  small  hydrometeors. 
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Hydrometeor  mass  loss  results  for  free  stream  Mach 
number  = 20. 


For  larger  shock  standoff  distance,  such  as  might  apply  to  off-stagnation  point  locations, 
mass  loss  of  the  particle  becomes  more  significant. 

Hydrometeor  slowdown  results  for  the  Mach  20  conditions  are  shown  in  Figure  2-11. 
The  figure  indicates  that,  although  mass  loss  is  small,  particle  velocity  may  be  changed 
significantly,  particularly  for  smaller  particles. 
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Figure  2-11.  Hydrometeor  slowdown  results  for  free  stream  Mach 
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SECTION  3 


CONCLUSIONS  AND  RECOMMENDATION 

On  the  basis  of  the  present  study,  the  following  conclusions  were  reached  relative  to 
the  prediction  of  mass  losses  from  droplets  immersed  in  reentry  vehicle  shock  layers: 

• Smaller  particles,  Rp  < 20  pm,  break  up  and  need  not  be  considered  in  a thermal 
response  analysis. 

• The  viscous  density  ratio  determines  the  characteristics  of  the  mixing  region  about 
the  droplet;  current  analysis  indicates  that  the  air  motion  can  be  decoupled  from 
the  water  motion. 

• For  flight  regimes  and  particle  sizes  of  interest,  the  se.ni-infinite  slab  assump- 
tion is  adequate  for  predicting  transient  surface  temperatures. 

• Transient  response  of  hydrometeors  need  not  be  considered  in  droplet  mass  loss  cal- 
culations unless  the  freestream  Mach  number  is  below  10. 

§ Heating  events  for  the  unblown  air  boundary  layer  about  the  droplet  were  obtained 
and  correlated. 

• Ablative  mass  loss  of  a hydrometeor  was  found  to  be  less  than  15  percent  of  the 
initial  mass  of  the  droplet  for  all  particle  sizes  and  Mach  number  conditions  of 
interest. 

In  sunmary,  these  analyses  have  indicated  that  rain  particle  with  radii  greater  than  20  urn 
remain  essentially  intact  a?  they  traverse  the  bow  shock  layer  in  the  stagnation  region  of  a 
high  speed  reentry  vehicle.  Particles  smaller  than  20  pm  radius  are  likely  to  break  up  from 
the  aerodynamic  pressure  in  the  shock  layer,  but  particles  larger  than  20  pm  will  survive 
either  because  transient  thermal  response  keeps  the  drop  below  vaporization  temperatures  (low 
Mach  number),  or  because  the  hydrometeors  do  not  lose  significant  mass.  Particles  in  the  20 
to  100  pm  radius  range  do  experience  some  slowdown  (relative  to  the  vehicle). 

Based  on  these  analyses,  it  is,  therefore,  recommended  that  coupled  hydrometeor  slow- 
down and  ablation  calculations  including  blowing  effects  be  performed  at  conditions  where 
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droplet  vaporization  Is  anticipated.  Experimental  verification  of  droplet  transient  response 
is  reconnended  for  particle  diameters  in  the  100  to  300  pm  range  at  Mach  number  greater  than 
10.  Such  testing  would  also  Indicate  the  effects  of  the  super-critical  surface  state  and  drop- 
let shape  distortion  effects  on  the  conclusions  reached  here. 
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APPENDIX  A 


ENVIRONMENTAL  TRAJECTORY  CHARTS 
by 

M.  D.  Jackson 

Reentry  parameters  for  a vehicle  a. id  a particle  entering  the  vehicle  shock  layer  are 
presented  in  this  Appendix.  The  parameters  are  defined  as  a function  of  velocity  and  altitude. 
The  range  of  velocities  and  altitudes  considered  was  30  to  0 kft/sec  and  60  to  1 kft. 

The  AC"  code  (Reference  A-l)  was  used  to  characterize  the  environmental  parameters. 
Calculations  were  performed  for  air,  with  ionization  and  dissociation  considered.  Table  A-l 
lists  the  species  which  were  used  in  the  calculations.  Diffusion  factors  (F.^)  were  also  used 
in  tho  computations,  and  are  defined  at  the  end  of  Table  A-l.  Freestream  environmental  pro- 
perties rf  the  atr  ' namely  density,  pressure,  temperature  ana  sonic  velocity  are  sum- 

marized in  Fig  Twi  f-l  A*2,  A-3,  and  A-4  respectively. 

Table  A-F  su.  .ar'.es  the  parameters  calculated,  their  associated  symbol,  units,  defin- 
ition and  correpondir  .ire  number.  The  following  parameters  listed  in  Table  A-2  are  re- 
lated the  environment  of  a particle  entering  the  vehicle  shock  layer: 

• particle  stagnation  pressure  and  total  enthalpy 

• unit  Weber  number 

• unit  Reynolds  number 

• particle  Mach  number 

In  the  ACE  code  it  is  assumed  that  gas  mixture  can  be  treated  as  a collection  of  ther- 
mally perfect  species.  Real  behavior  of  gas  mixtures  become  important  as  the  mixture  approaches 
its  critical  state.  However,  for  the  pressures  and  temperatures  presented  herein,  these  ef- 
fects are  insignificant. 
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Figure  A-2.  Freestream  pressure 
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Freestream  temperature. 
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Figure  A-4.  Freestream  sonic  velocity. 


TABLE  A-2.  SUfYIARY  OF  TRAJECTORY  PARAMETERS 


ooo* oe 


CM 


(D0S/*aaj)  ApotaA 


A-8 


Figure  A-5.  Freestream  Mach  number  (M, 
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Figure  A-13.  Particle  stagnation  pressure. 
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APPENDIX  B 


DETAILS  OF  LOCK'S  FORMULATION 

On  the  usual  assumption  that  the  change  of  velocity  from  U2  to  takes  place  in  a 
layer  of  small  thickness,  and  that  v is  everywhere  small  compared  with  u,  the  boundary  layer 
equations  are 


u 


+ v = 
3x  9y 


for  the  upper  fluid,  and 


u 


1“  + v- 
ax  v ay 


au  . 


(B-l ) 


(B-2) 


for  the  lower  fluid,  where  and  v2  are  the  kinematic  viscosities  of  the  two  fluids. 
The  equation  of  continuity  is 


3u  + 3v 

ax  ay 


0, 


so  that  there  exists  a stream  function  ^ such  that 


u 


v 


ax  * 


In  order  to  solve  Equation  ( A- 1 ) we  use  the  nondimensional  variable 


and  look  for  a solution  in  which 

1>  = (v-jU-jX) 


(B-3) 


(B-4) 


(B-5) 


B-l 
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Then 


u = U-jf | (ri| ) . 


^l^r)k  ^(n,)  - f,<V» 


(B-6) 


(B-7) 


and 


U, 


l/2 


ly=  U1  ( v^x)  fl(nl}  * 


(B-8) 


Equation  (B-l)  then  reduces  to 


d'f,  d2f 

2 W f,  L = 0 

dn|  dn| 


(B-9) 


In  the  lower  fluid  it  is  convenient  to  use  a different  variable  n2  defined  by 

- ( U]  t y 
n2  ‘ \^y)  ' ’ 


(B-10) 


so  that 


n2  = \ V * 


and  to  put 


\\>  = (vgU-jX)  'z  t2(n2) 


(B-l 1 ) 


Then 


u = U1f2(n2)  » 


/ \ ‘/2 

v * f(nr)  ,n2f2(n2)  ' f2(n2)) 


(B-l  2) 


(B-13) 


and 


a».„  ZL 


. V: 


§ = MvJ  w- 


(B- 14) 
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Equation  (B-2)  then  reduces  to 


d3f,  d2f, 

2 — £ + f — L = o . 
dn^  1 dn£ 


(B-15) 


Boundary  Conditions 

The  boundary  conditions  at  infinity  are 


u ■*  as  n-j  -*■  +0° 


and 


u -►  as  n2  -00 


so  that 


fj  -*•  1 as  ri|++® 


(B-16) 


and 


f2  * 07  35  n2 


(B-17) 


Since  the  motion  is  steady,  the  interface  between  the  fluids  is  streamline  ->  = 0 which 
passes  through  the  origin,  and  is  therefore  given  by 

V f2  = 0 • 

If  is  the  value  of  such  th^t  f^njj)  = 0,  then  we  must  have  also 


f^njj)  = 0 . 


where 


0 _ v2 


n2  = l v1 


'2 


The  other  boundary  conditions  to  be  satisfied  at  the  interface  are  that  the  velocity 
and  the  normal  and  transverse  components  of  stress  should  be  continuous. 
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APPENDIX  C 


HYDROMETEOR  ACCELERATION  AND  MASS  LOSS  IN  A SHOCK  LAYER 
CONSIDERING  BLOWING  AND  REAL  GAS  EFFECTS  ON 
HEAT  AND  MASS  TRANSFER 


Jaffe  (References  C-l  and  C-2)  coupled  the  mass,  momentum,  and  energy  relations  for  an 
ablating  particle  in  the  shock  layer  of  a high  speed  vehicle.  This  development  is  here  mod- 
ified to  include  the  effects  of  blowing  and  real  gas  effects  on  droplet  heat  and  mass  transfer 
In  addition,  a more  appropriate  dependence  of  heat  transfer  coefficient  on  shock  and  particle 
conditions  is  employed. 


Particle  Momentum 

Force  on  a particle  in  the  shock  layer  of  a high  speed  vehicle  is  given  adequately  by 
a drag  coefficient  formulation  as  follows: 
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where 

CQ  = drag  coefficient* 

Ac  * particle  cross  sectional  area,  ft2 
P2  = shock  gas  density,  lbm/ft3 
u = particle  velocity,  ft/sec 

For  a spherical  particle 
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For  a sphere  at  hypersonic  conditions,  Cq  = 1.0  approximately. 
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Hydrometeor  Ablation 

Ablation  of  a hydrometeor -particle  can  be  evaluated  by  considering  mass  and  energy  con- 
servation in  combination  with  boundary  layer  and  droplet  properties.  Mass  and  energy  control 
volumes  are  illustrated  in  the  following  sketch. 
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Mass  conservation  gives 
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where 

ti  = average  surface  mass  loss  flux 
Ap  = particle  surface  area  = 4nRp2 

The  average  ablation  flux,  m,  can  be  expressed  in  terms  of  the  stagnation  point  values  under 
the  following  conditions 


C-2 


• The  leeside  hemisphere  experiences  low  heating  and  does  not  lose  mass. 

• The  average  windside  hemisphere  heat  transfer  and  ablation  rate  is  71  percent  of 
the  stagnation  point  values  (see  Section  2). 

That  is, 

m = 0.71  m0jSp/2  (C-6) 

The  quantity  m^  <.p  is  the  stagnation  point  mass  removal  flux  and  can  be  evaluated  from 
energy  conservation  at  the  particle  stagnation  point. 

In  order  to  evalaute  m^  <-p,  the  control  volume  is  chosen  to  be  a thin  layer  extending 
from  slightly  above  the  ablating  surface  to  beneath  the  thermal  penetration  depth.  This 
energy  control  volume  moves  with  the  receding  surface  at  a rate  equal  to  the  surface  recession 
rate.  Energy  conservation  for  this  control  volume  gives  the  following 

Q --  VHW  - Hf)  ♦ f = PeueCH(He  - Hw)  (C-7) 

where 

dE/dt  = the  rate  of  energy  storage  = 0 in  steady  state 
n^C^  = real  gas  film  coefficient  including  blowing  effects 
Hg,  Hw,  Hf  are  enthalpies  at  the  boundary  layer  edge,  the  surface  of  the  droplet,  and 
the  initial  temperature  of  the  droplet,  respectively 

Equation  (C-7)  is  useful  when  ablation  and  heat  conduction  reach  steady  state  conditions  and 

when  the  thermal  penetration  depth  is  much  less  than  the  droplet  radius.  A crude  estimate  of 

thermal  penetration  before  vaporization  begins  is  given  by 

-d  = -VTVAR-Ti)/%  <C"8> 

For  Rp  = 20  um,  Mw  = 20,  z = 15  kft 

d = (0.88  x 1.0"4  Btu/ft-sec°R)(1165  - 400)/(10  x 106  Btu/ft2sec) 
d = 0.67  x 10'8  ft  = 2 x 10'3  pm 

Thus,  only  a very  thin  surface  layer  of  the  20  pm  hydrometeor  is  affected  by  the  heating  prior 

to  vaporization.  Similar  arguments  show  that  thermal  penetration  and  storage  rate  are  neglig- 
ible when  vaporization  occurs.  Equation  (C-7),  thereto  reduces  to 
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3>  _ <He  - Hw> 

peueCH  <Hw  - Hf> 

Equation  (C-9)  can  be  combined  with  the  blowing  effects  correlation  described  in  Section  2 to 
give  itiq  in  terms  of  enthalpies,  the  non-blown  transfer  coefficient  and  the  blowing  parameter, 
X. 

That  is,  from  Dorrance  (Reference  C-4), 
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where  subscript  "o"  denotes  non-blown  values.  Equation  (C-10)  can  be  rearranged  to  give 

peuecH 

mD  s 2x  ^n(l  + 2Am0/peueC^)  (C-ll) 

Combining  Equation  (C-9)  with  Equation  (C-U)  and  applying  the  result  at  the  droplet  stagna- 
tion point  yields  the  desired  result 
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At  the  hypersonic  conditions  of  interest,  the  proper  boundary  layer  edge  enthalpy  (Hg)  for 
use  with  a stagnation  point  film  coefficient  is  the  total  enthalpy  behind  the  second  shock, 

Htot*  this  enthalPy  is  generally  large  compared  to  the  wall  enthalpy,  Hw<  Also,  the  enthalpy 

difference,  Hw  - Hf,  is  the  energy  absorbed  in  heating  and  vaporizing  a unit  mass  of  the 
hydrometeor.  The  quantity  has  been  termed  the  "heat  of  ablation",  Qp,  and  would  include  a 
heat  of  fusion  if  the  hydrometeor  were  initially  an  ice  particle.  Approximations  for  the 
enthalpy  terms  can  be  used  because  the  mass  flux  varies  as  the  logarithm.  For  convenience. 
Equation  (C-12)  is  expressed  as 

peueCH  r , 
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Combining  Equations  (C-5),  (C-6),  and  (C-13)  and  substituting  the  mass  of  the  particle  gives 


C-4 


(C-14) 


dM. 


dRn 


un  j 

3t^  * 3t  (^RrVd)  * 4nRpJPP  af 

°-7,peueCHc]5p  l"  ['  + m'toM 


* -(4ttRp2) 


AX 


or 


dRp  °'71peueCH  ]sp  *"  f1  + Kot'Op] 

ar <c-15> 


To  this  point  the  development  of  the  coupled  mass,  momentum,  and  energy  equations  de- 
scribing the  response  of  the  hydrometeor  is,  except  for  the  addition  of  blowing  effects,  quite 
parallel  to  that  of  Jaffe  (References  C-l  and  02).  Contrary  to  Jaffe,  it  is  felt  that  for 
particles  in  supersonic  flow,  the  heat  transfer  coefficient  should  be  scaled  as  the  particle 
velocity  divided  by  the  square  root  of  the  particle  radius.  This  has  justification  in  the  theory 
of  Fay  and  Riddell  (Reference  03)  and  in  the  numerical  work  performed  by  Bartlett  and  Putz  (Re- 
ference 04).  Thus, 
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Jaffe  (References  C-l  and  02)  utilized  a Nusselt  number,  based  on  initial  conditions,  which 
varied  inversely  with  particle  radius  and  independent  of  velocity.  BLIMP  calculations  and  the 
correlation  given  in  Section  2 provide  the  initial  condition  heat  transfer  coefficient, 
peueCHo  SP  i*  Equation  can*  therefore,  be  written: 
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Normalizing  Equation  (017)  and  converting  to  path  length  as  the  independent  variable  gives 
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Boundary  condition  on  Equation  (C-17)  is  R = 1 at  x * 0.  Integration  gives 
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For  the  spherical  particle 
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At  the  nosetip  surface  where  x = 1.0 

H^T  = (’  ' I Kf)  (C‘21) 

In  Equation  (C-2I)  must  be  less  than  2/3  for  the  particle  survival  to  the  nosetip 
surface.  greater  than  2/3  corresponds  to  complete  particle  ablation  in  the  shock  layer. 
Equation  (C-18)  shows  that  the  particle  radius  reduction  is  not  a significant  function  of  the 
particle  slowdown  since  the  drag  constant,  Kp,  does  not  enter  the  relation.  Some  elements  of 
the  energy  constant,  K£,  such  as  H^ot,  are  dependent  on  the  gas/particle  relative  velocities. 
Thus,  for  conditions  of  significant  slowdown,  Equation  (C-21),  for  example,  would  become  in- 
accurate. The  equations  are  accurate  for  conditions  of  interest,  however. 


Coupled  Ablation  and  Slowdown 

Combining  Equation  (C-19)  with  (C-4)  using  the  boundary  condition  u 8 1 at  x = 0 yields 


u 8 exp 


(C-22) 


where  > 2/3  for  particle  demise  in  shock  layer.  Parametric  results  using  Equations  (C-19) 
and  (C-21)  are  shown  in  Figure  C-l. 

Hydrometeor  mass  and  velocity  attenuation  by  the  technique  presented  above  is  somewhat 
less  than  the  results  obtained  using  the  technique  of  Oaffe  (Reference  C-2),  which  gives 
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Results  are  not  drastically  different  because  Nusselt  number  correlations  give  less  heat  trans- 
fer than  predicted  for  real  gas  conditions  by  BLIMP  whereas  the  BLIMP  heat  transfer  results 
have  been  reduced  for  blowing  effects. 
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